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I. INTRODUCTION 


Since Kolmogorov derived asymptotic results for the second order structure function 
(Am^) in the dissipation range and for the third order structure function (Au^) in the iner¬ 
tial range in 1941 from the Navier-Stokes equations under the assumption of homogeneity, 
incompressibility and isotropy, a vast amount of work was carried out on the subject. One 
of the reasons is that the scale of with the viscosity u and the mean dissipation (e) in these 
results was thought to be universal for all kind of flows. Soon, it has been found that higher 
orders do not follow the scaling predicted by Kolmogorov and the deviations from it be¬ 
come stronger for higher orders. Hill [l| and Yakhot independently derived higher order 
struture function equations. However, the source terms are unclosed. For more insight, 
we derive equations for the source terms. Applications will be presented in papers to be 
submitted and the present article is meant as a basis for these. 


II. FOURTH ORDER STRUCTURE FUNCTION EQUATIONS 


We start be deriving the equations for hi^er order structure functions, the derivation 
of which has already been carried out by Hilljj, |^. The reason for re-deriving them here is 
twofold. Firstly, the present article serves as basis for papers to appear in the refereed jour¬ 
nals and for that purpose it is reasonable to have all equations in one document. Secondly, 
additional steps used in the derivation of the structure function equations are required for 
the derivation of their source term equations. 

Only the steady state form of the averaged equations will be considered. The two points 
are denoted by a; = {xi,X 2 ,X 3 ) and x' = {x[,X 2 ,x'^). Assuming incompressible flow, the 
momentum equations in the Navier-Stokes equations are written at the two points 
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dt dXn 
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Here Ui and u' are the components of the velocity, p is the pressure and p the kinematic 
viscosity. Einstein’s summation convention for indices appearing twice is used. The density 
was set equal to unity. These equations are completed by the continuity equation which 
holds at both points. 
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Substracting eq. ([2]) from eq. ([T]) one obtains for the velocity increment Aui defined by 
Aui = Ui{x) — u'^{x') the equation 
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Here the difference of the pressure gradient at the two points is dehned as APj and the 
pressure is given by the Poisson equation 


d'^p dui duj 
^ dxj dxi 


= 0 


( 6 ) 


From this an equation for the pressure gradient dp/dxk may be derived by diferentiation 

f dp 


dxi \dxi 
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Rk 


which defines the operator Rk. The equation for AP^, multiplied by the viscosity, then 
becomes 

/a^APfe , a^APfe 
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-^ + z/ (Pfc — P^) — 0 


V dxi ' y ' ^ ^ ^ 

Equation (j5]) is then multiplied by Awj and the corresponding equation for Awj is multiplied 
by Awj to obtain 
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Adding eq. ([9]) and eq. flTOj) one obtains an equation for AuiAuj, 
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This also defines the instantaneous values of the dissipation Sij at x and eh at x'. Equa¬ 
tion fITT]) is an equation for velocity differences at two points. However, the derivatives do 
not reflect this and are still carried out with respect to x and x'. Changing the independent 
variables from x and x' to the new independent variables 


X = -{x + x'), r = {x - x') 


and using the transformation rules 
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(15) 


one obtains for the transport term in eq. mil 

dAuiAuj , dAuiAuj dAunAuiAuj 


Ur. 


dxr 


+ Ur 


dxL 


dvr 


where the derivatives d/dXn were neglected, as they vanish after averaging due to the 
assumption of homogeneity. Likewise the Laplacian in eq. ffTTll ) becomes 


d^AuiAuj ^ d'^AuiAuj ^d^AuiAuj 


+ 


(16) 


dx"^ dx'^ dr^ 

We now average eq. mil. Dehning the second order longitudinal structure functions as 

^ 2,0 = (Aul) (17) 

and the transverse structure function as 


^ 0,2 = (Aul) (18) 

and by setting i = 1, j = 1 and i = 2 , j = 2 , respectively, in eq. m» one obtains the 
transport equations for 5*2,0 and 5*0,2• We can further simplify these equations under the 
assumption of isotropy. Specihcally, we can hnd the isotropic form of both the gradient 
and the Laplacian, so that they depend on the separation distance r only, instead of the 
separation vector r„. For the steady state case after transforming the transport term and 
the Laplacian into their isotropic form (cf. Robertson jsl), we have 


<95^3,0 I 2 ^ 4 _ 
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or r J r 

5 * 5*1 2 ,2 I 

-h - 5 i ,2 I = 2 i/ 

or r 


5^ .52,0 2 I ^ 

r., n -1-*J2,0 H- ^ Wo,2 — *->2,0j 


d^Si 


0,2 


dr'^ 


H—*5*0,2-5- (*5*0,2 — *5*2,0) 


These are the second order structure function equations given by Hill 
noted that due to isotropy 

(£11) = (£22) = 3 • 

Continuity provides an additional relation between 5 * 2,0 and 5*o,2 

r 55*2,0 


2 dr 


+ *5*2,0 ~ *5*0,2 — 0 


as well as between S 3 0 and 5*i 2 


-2{eu) (19) 

- 2 (£ 22 ). (20) 

j^. In Q it is also 

( 21 ) 

( 22 ) 


r + 5*3 0 — 65*1 2 — 0 (23) 

or 

Using these two equations in eq. flT^ . this leads then to the well-known Kolmogorov equa¬ 
tion IgI, 



(75*2,0 

dr 


4(e). 


(24) 


Kolmogorov had derived two asymptotic solutions of this equation, namely 5*2 ,0 = 
(e)r^/(15z/) for r —)■ 0 and 5 * 3,0 = 4/5 (e)r in the inertial range. The latter is known as 
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the four-fifth law. Next we consider equations for structure functions and equations for 
source terms at the fourth order. The procedure is the same as above, i.e. we multiply the 
equation for Aui (eq. (j5])) by AujAukAui and add up all four combinations resulting in 
an equation for AuiAujAukAui. In the next step, the coordinate transform outlined above 
is used and the gradient and Laplacian are re-written in their isotropic form. This is pre¬ 
sented in Appendix A below in some detail. Hillfl gave the necessary coefficients up to the 
seventh order and presented a matrix algorithm which can be used to derive the isotropic 
form for all orders. Note that the coefficients in the Laplacian of the higher order structure 
function equations given by Hill and consequently the matrix algorithm were corrected in 
http://arxiv.org/abs/physics/0102055, In their isotropic form the equations then read 
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(25) 


(26) 


(27) 


Note that we have three coupled equations for the three unknowns S' 5 , 0 , -S' 3,2 and S' 1 , 4 . 
Noticeably, this is the case for all even orders, but not the odd orders. The instantaneous 
values of the source terms are given by 


T 4 ,o = AAu\APi (28) 

T' 4,0 = 6An^ (£11 -|- (29) 

T2,2 = 2 AM 2 AM 1 AP 1 -|- 2Au\Au2AP2 (30) 

T' 2,2 = (^11 + ^n) + 4 A'UiA'U 2 (^12 + ^ 12 ) + (^22 + ^ 22 ) (31) 

To,4 = AAulAP2 (32) 

Eo,a — QAu\ {s22 + ^ 22 ) • (33) 


Different to the second order equations, there are now pressure source terms Tj,j ~ 
Ati*+i-iAP acting on the system of equations. They appear at all orders higher than 
the second, where they had dropped out due to isotropy (cf. Hill [^). Additionally, there 
are the dissipation source terms Ei,j. Note that again the second order is special, inasmuch 
as the averaged dissipation source terms equal the one-point quantities {si,j). For higher or¬ 
ders, this is not the case, i.e. Pj,j is now clearly r-dependent. Consequently, it is not possible 
to find asymptotic solutions in the spirit of Kolmogorov without closure assumptions. 

Therefore, it is worthwhile to study the higher order source terms in more detail, as 
anomalous scaling of the structure functions most likely originates from the source terms. 
For that reason, we proceed to derive the equations for all six source terms in the following. 
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III. FOURTH ORDER PRESSURE SOURCE TERMS 


Equations for the pressure source terms of the fourth order structure function equations, 
which have the generic forms UiUjUk^Pi, will be derived next. For that purpose we hrst 
derive an equation for the unaveraged product of three velocity increments by extending the 
procedure which led to eq. ffTTD to the product of three velocity increments. The result is 
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We multiply this equation by APi and convert the transport term as 

. „ dAuiAujAuk , . „ dAuiAujAuk dAuiAujAukAPi , dAuiAu^AukAPi 
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Ur,.— - \- U ' 


dXr 


dxi 
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(35) 

where the second term is interpreted as a sonrce term and is shifted to the right hand side 
of the eqnation. Adding eq. ([7]) mnltiplied by AuiAujAuk the eqnation then reads 

, dAuiAuj AukAPi 
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- AuiAPi [sjk + e'jk) - AujAPi {sik + e[f,) - AukAPi [sij + eh) 

The nnsteady term has been retained and was shifted to the right hand side of this eqnation 
because it does not vanish in the steady state case after averaging. The transport terms and 
the diffnsion terms in this eqnation are transformed into their isotropic form. For the trans¬ 
port term, we can not simply use the coefficients given by Hill, dne to different symmetries of 
the tensor AukR^ compared to Sij. This derivation is somewhat laborious and reported in 
Appendix B. Conseqnently, there are more scalar fnnctions needed to determine the trans¬ 
port term. Specihcally, two additional terms (Am 2 T 3 ^i) and (Am 2 T 3 _i) appear. These are 
dehned as 

(Am2T3,i) = {Au2 [3AuIAu2APi + AUIAP 2 )) (37) 
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and 


(38) 


(Am 27 i, 3 ) = (Am2 (AM1AM2AP2 + 3AM2AP1)) . 
(AmiT 4 ^o) is given by multiplying eq. fl^ by Ami and averaging, i.e. 


(AuiPj^o) = 4 Am^APi 


(39) 


and similarly for (AmiT 2 ^ 2 ) and (AmiTo, 4 ). The Laplacian however has the same form as 
given by Hill, so the coefficients can be carried over. By setting the indices according to the 
dehnitions in eq. (125]) . eq. (150]) and eq. (152]) we obtain after averaging the following equations 
for ( 24 , 0 ), ( 22 , 2 ) and (To^ 4 ) in the isotropic form 
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d (Atiir„,4) I ^ i (AujTi.a) = - (12At.2A«2AP2Ai’2> 
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IV. FOURTH ORDER DISSIPATION SOURCE TERMS 


The dissipation source terms in eq. ([29]), eq. (l3l]) and eq. fl33|) have the generic forms 
AuiAuj{eki + e'f^i). In order to derive the equations for these terms we first derive an equation 
for {ski P^'ki)- An equation for the instantaneous value of the product of velocity derivatives 
is derived by first taking the derivative of the momentum equation written for Uk with 
respect to Xm and combining it with the corresponding equation for the derivative dui/dxm 
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Noticeably, 2i'Xij can be interpreted as the dissipation of Eij. Using the definitions of Eij 
and in eq. (ITT]) we obtain 
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dun duk dui ^ du'^ du'f. du[ 
dXm. dXrt. dXrn. 


dx'^ dx'^ dx'^ 




2 z/ ( 


dXm dXn dXr. 


jm 

V 

J^lk+A', 


dx'm dx'^ dx'^ 


Ik 


2 ^ duk d'^p ^ du'j^ d'^p' 


dxm dxidxm dx'^ dx[dx'^ 


Pkl+Pki 

^ I dui d‘^p ^ du'i d'^p' 


dxm dxkdxm dx'^ dx'^dx'^ 


Plk+Plk 

^ I ^^gfcZ + ^'kl 


dxl 


n 


dx'^ 


2v {Xki + xlz) 


(46) 


which also defines the quantities A^i + A'^^i and P^i + P^i- Combining this with the equa¬ 
tion ffTTD one obtains 
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(47) 


The term (Eij + E(j)(Eki + e'^^i) in this equation will, after averaging and using the definitions 
in eq. fl29|) . eq. fl5T|) and eq. fl33|) generate the new dissipation parameters at the fourth 
order (£ 12 ), 411 ^ 22 ) and {E 22 ) discussed in the main text. The isotropic form of the 

transport terms in the equations for ( 772 , 2 ) and (i 7 o, 4 ) are derived in Appendix B 
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using the method outlined by Robertson [^. They are 
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5(A-Ui£'o,4) ,2/a pi \ \ — 
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- (^ 0 , 4 ) 

— 24z/ {A.U 2 ^U 2 (-P 22 + P 22 )) 
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- 12 {uAu 2 Au 2 (X 22 + X 22 )) 

- 6 {^(£22 + ^ 2 )^) 
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(55) 


V. SIXTH ORDER 


We also present the equations for the sixth order. The reason for doing this is that the 
source term of the dissipative source term of the sixth order contains the triple product 
{eijeklsmn) 1 he. we hud the third moment of the dissipation in the system of equations. In 
fact, all moments of the dissipation are found in the system of equations, when one continues 
to derive the source terms at higher orders. 

At the sixth order the structure function read according to Hill (http: //arxiv. org/abs/physics/0102C 
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(I- +-) ^56 = - (Ae) - (Bo,6) 

\dr r J 
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Tg^o =6 Ami^APi (60) 

Eefi =15Ami^ (^ii + (61) 

T4,2 =4Ami3Am2^APi + 2Ami^Am2AP2 (62) 

-£'4,2 =6A'Ui^A'U 2^ (sil + ^ii) + 8 A'Ui^A'U2 (^12 + ^ 12 ) 4“ Ami^ (^22 + ^ 22 ) (63) 

T2^4 =2Ami Am2^APi + 4Ami^Am2^AP2 (64) 

p2,4 =A'U2^ (^ll + Ell) 4“ SAttiA'U2^ (^12 + ^ 12 ) 4" 6A'Ui^A'U2^ (^22 + ^ 22 ) (65) 

To ,6 = 6 Am 2 ^AP 2 ( 66 ) 

Eofi =15Am2^ (^22 + ^ 22 ) (67) 


We will only consider the dissipation source terms because those will generate the third 
moment of the dissipation distribution at the end. In order to derive equations for these 
source terms we hrst derive an generic equation for the fourth product of instantaneous 
velocity increments to be called 


Au"^ = AuiAujAukAui 


( 68 ) 


Combining eq. (A. 12) for AuiAuj with its form for AukAui one obtains 


dAu^ 

dt 


dAu^ 


Ur 


OAu^ 

u'„ ^ = —AujAukAuiAPi — AuiAukAuiAPj 


dxn "■ dx'^ 

— AuiAujAuiAPk — AuiAujAukAPi 
fd^Au^ d^Au^ 


- 2z/ 


dx'^ 

(dAuiAuj dAukAui dAuiAuj dAukAui 


(69) 


\ dXn dXr. 
AukAui{£ij + eh) 
AuiAuj{eki + e'ki) 


+ 


dxL 


dxi 


This will be combined with eq. fl46|) to obtain for the combination of terms in eq. fl6Tl) . 
eq. fl6^ . eq. fl6^ and eq. fl^ for the generic form 


Au^ [Spq + Epq) = AuiAujAukAui [Spq + e'pq^ 


(70) 
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the equation 


<9Am^ {spg + e'pg) 


dt 


^ ^ dAu^epg + e'pg) ^ ^ , + 

~r Un “T 


dXr. 


dxL 


— {AujAukAuiAPi + AuiAukAuiAPj + AuiAujAuiAPk 
+ AuiAujAukAPi) (spq + e'pg) 

— 2uAu^ (^pg + ^'pq + ^gp + ^gp) 

— 2uAu {Ppq + Ppq + Pqp + Pqp) 

( d'^Au^ [Spq + e'pq) d'^Au^ [Spq + e'pq) 


- 2v 


\ 9x2 

^Au"^ de 


pq ' ^pq 


dXn dXn 

- 2uAu^ {Xpq + x'pq) 

— 2v {e. 


+ 


dx'^ 

dAu^ de. 


4 - 

pq ^ ^pq 


dx' dx' 


'pg 


, ^ f dAuiAuj dAukAui dAuiAu^ dAu^Aui 

+ Sg)(^^--+ 


dxL 


dxL 


(71) 


dXn dx. 

— AUkAui {sij + e^^ (^Epq + Spq^ 

— AUiAUj {Skl + £kl) ippq P ^pq'} 


From this form one could derive equations for (Eq^), (-^ 4 , 2 ), (-£ 12 , 4 ) and (Eq q). We will not 
write those equations down but want to consider a generic equation for the last term on the 
right hand side of eq. (1711) to be called 


Eijklpq 15 AUj Alij (Sfc; -|“ ^fc;)(^pg P ^pq')' 


(72) 


For this purpose we hrst derive an equation for = {£kiPPki)i^pgP^pg) from eq. which 


reads 


de^ 


P Ur 


9^2 


dt ' "dXr 


, 9g2 
dxL 




- 2z/ {spq + e'pq^ {Aki + + Aik + 4^^) 

- 2z/ {Spq + e'pq^ {Pkl P Pkl p Plk p Plk) 

- 2l^ + e'ki) {Apq + A'pq + Aqp + Agp) 

- 2l^ {Skl P e'kl) {Ppq + Ppq p Pqp + Pqp) 

de^ de^ 

P 


2u 


9x2 Q^/2 

dski P e'p.1 dspq + dski P e'ki dspq + e,. 


'pq 


+ 


-pg 


dxL 


dxL 


(73) 


dXn dXn 

— 2z/ {Skl P e'ki) {Xpq P Xpq) 

- 2z/ {Spq + e'pq) iXkl p x'kl) 

Combining this equation with eq. flTT]) one obtains after averaging for {AuiAuje"^) = 
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{AuiAuj{eki + £'ki){,£pq + e'pq)) the form 

, dAuiAujS^ 


dAuiAunE^ 

dx. + 


dxL 


- {e^ {AujAPi + AuiAPj)) 

- 2v (^AuiAuj [spq + e'pg) {A^i + A[i + Aik + ^zfc)) 

- (^AuiAuj {spq + e'p^ {Pki + P'ki + Pik + PL)) 

- 2v (^AuiAuj {ski + sL) {Apq + A'pq + Aqp + ) 

- 2v l^AUiAUj (Skl + e'kl) {Ppq + Ppg + Pqp + Pqp) ) 

/ / AuiAujS^ AuiAujE^L \ 

2 / f dAuiAuj de'^ dAuiAuj de^ \ \ 

^ \ V dxn dxn ^ dx'^ dx'^) / 

dXn dXn dx'^ dx'^ 

- 2v {AuiAuj {ski + s'ki) {xvq + XpJ > 

- 2i/ {AuiAuj {spq + ixki + x'ki)) 

- ([eij + e'ij) {ski + sL) {^pq + Spq) ) 



2u (^AuiAuj 


(74) 


In this equation the term + e[j){eki + £$,;)(£pq + 6^^)) will generate third moments of 
the dissipation distribution. The longitudinal form {Fe^) = 15 (Awf (en + e'n)"^) of eq. fTMll 
will contain the term 15 {(en + ^ii)^) which generate parameter for large r. 


VI. DISCUSSION 

The procedure of deriving equations which parameters proportional to (e") is based 
on identifying source terms where derivatives square are multiplied by AuiAuj. A hrst 
example is the fourth order dissipation source term = 6Aul{eii + e'n) in eq. (12^ 
which led to the term 6 ((^n T^n)^) in eq. (USD and hnally to the dissipation parameter 
Other source terms which could be considered are AuiAujAPkAPi in eq. fl5B]) or 
the terms 2uAuiAuj{Aki + and 2i'AuiAuj{xki + x'ki) i'^ In the equation for 

AuiAujAPkAPi parameters {sij + e[j)APkAPi then appear which are proportional to veloc¬ 
ity derivatives squared times pressure gradients squared. In the equations for the F-terms 
(eq. d72|) I we hnd terms which are the product of hve velocity derivatives. Those become 
constant in the inertial range but not in the dissipative range. Their contribution to inter- 
mittency are not clear at this stage. 
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Appendix A: Laplacian of fourth order structure functions 

Following Robertson j^, a fourth order tensor of two-point type which is invariant to 
rotation and reflection of the coordinate system (cf. isotropic turbulence) is given by 


^ Tin , ^ ^ r,rk , , ^ TiTj 


TirnThTl ruT] .,,,,, , 

^ijkl = -7-1“ ^2 %—2 ^ —7 1“ - 5 h A^6jk —^ + A^djl —5 h AfSkr 

+ Agdijdki + AgdikSji + AigdjkSii, 

(Al) 

where A* are scalar functions, 6 ij is the Kronecker delta, i.e. 6 ij = 1 for i = j and 6 ij = 0 
for i j and r* a separation vector with magnitude r, i.e. Xi = x[-\- r* where Xi and x' are 
two points in space. 

In the following, let Aijki = Sijki = {AuiAujAukAui), where Au = Ui — u[, i.e the fourth 
order structure function. As Sijki = Sjiki = Sijik = Sjuk = ..., A 2 = A^, = ... = Ay = S 2 and 
As = Aq = Aio = S 3 , i.e. from eq. flAip 


riTjrkVi 


TkTl 


Sijki — Si —-h S2 ( Sii ——h 6 ik^-Tr A Sii -^—|- Sjk —7- J- S 


rjVk 


TiTi 








TAk , . TiVj 

+ Okl^ 


A S 3 (^Sijdki A SikSji A SjkSiij. 


(A 2 ) 


Next, the scalar functions Si, S 2 , S 3 need to be determined. Without loss of generality. 


let ri = r, r 2 = rs = 0 . 

Choosing Ami, eq. fjA2p then yields 



Ami = A4 ^o = Al J- 6 A 2 A 3 A 3 , 

(A3) 

while A 1122 gives 


A 1122 — A 2,2 — A 2 J- As 

(A4) 

and A 2222 


A 2222 — Ao ,4 — 3As. 

(A5) 


Thus, the three scalar functions Si, S 2 and S 3 are determined by the three tensor components 
Aim, A 1122 and A 2222 and solving for them gives 


Al 

A 2 

As 


Aim — 6A1122 + A2222, 
A1122 — 2A2222, 

-A2222- 


(A 6 ) 
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Next, d'^Sijki/dr'^ depending on the scalar functions Si, S 2 , S 3 and the separation vector 
Ti is derived, for which the following relations are needed 


d 

dVn 

d 


Vr, 

y,2 

'r, 


dVr, 


r 

'r. 

r 


= 1 

= (Sin - 

1 


TiTr, 

^2 




2 

r 


dA{r) Tn dA{r) 


dr. 


dvn 
d (\ 

r 


r dr 
1 


dr. 


drr, 


- - — 

(diri. r» 


^2 ^ 


(A7) 


Consequently, 


Tn d 
r drr, 


( rr/.Si, 


rfY' . fY' 

n’ i’ n 


Using eq. flA7p and eq. flASp then results in 


1 / ,1 
-y = (r-i - rj -y = 0. 


(A 8 ) 


drl 


f c< ^i^fkri\ rirjrkri 

r' r^ )~ 


d^Si 2 dSi 


dr‘i 


;Si 


+ 


r 

Tjri 

.2 


rkTi 


i r^ 


Sjk 


Sij - 
Tjrk 


r^ 

8 „ rirjrkri 

—• 


r dr 

Tirj 


+ 




Tirk 


nrk' 


Sji - 


r,ri 


r.2 J 




+ 




Ski — 


rkri\ 


rfr\ 

r2 J (A9) 


Similarly with eq. flA7p and eq. flASp 

d^ r 


^ n 


So s, 


rkri , e- rjri r^rk 




+ s. 


ik ■ 


+ Sii- 


+ + Ski 




= S, 


fyi ^ 




+ 6,i^ + 6ki^ 


d^S2 2 dS2 

+ 


+ 2 h 

~ ^2 

+ Sii 
+ Sji 


Sij 


Sjk 


Ski — 
Tjrk 


dr‘i 

rjn\ _ 

r 2 ) 


r dr 
Jkri\ 


- 2 


Tjrk 


r2 J 

+ Sjk 


+ Sik 

Sii — 


Sji - 

Tiri 


r,ri 


- 2 


- 2 

rfi 


rjri 


S^k - - 2 ^) + Ski 


. nrj 

^2 


- 2 


nr. 


(AlO) 


Finally, 

d^ 


dr‘i 


d'^S 2 dS 

[S 3 {SijSki + SikSji + SjkSil)] = {SijSki + SikSji + SjkSil) I ^^2 + 


(All) 
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Now, adding eq. flA9p . eq. flAlOj) and eq. flAllj) and snbstitnting i = j = k = l = l gives 


drl 


Si 2 dSi 
- - H- - 

^ 'Y' 


:d‘^S 2 2882 

61 -:r^ + - 


dr'i 

d^S. 


+ 


r dr 
2 dSs 


-^Si 


24 

^*->2 


9r2 r dr 


Using the relations eq. flA 6 p yields 


d^Snn fd^Suu ^ 2dSini\ 8 ^ , 24 ^ 

-- — I -- 1 ---H-y*J 1122 

Yy'Y'^ 'Y' J ry'^ 


(A12) 


(A13) 


Setting i=j = l,/c = / = 2in the summation of eq. 


, eq. flAlOl) and eq. flAllI) yields 


d^S 


1122 


dr'i 

n 


2 ^ ( 8^82 28 S 2 


8'^ S3 ^ 2 8Ss 


2 o 

- n 02 


(A14) 


and with the relations eq. 


^51122 _ fd'^Si 


(A16) 


In the same way, setting i = j = k = l = 2 gives 

^'^82222 _ 12 10/ 2 883 \ 

( 9 r^ j ,2 2 ^ Q ^2 r 8 r ) ' 


i.e. using eq. flA 6 p 


8‘^S2222 / 8“^82222 2 8S2222 \ , 12 4 

- 7^ - = - 7^ - 1 -- H - 77O1122 - 5 -> J2222 

8 r^ \ 8 r^ r dr J r^ r^ 


(A16) 


(A17) 


These are the same equations Hill derived using a matrix algorithm (note there is a re¬ 
cent correction to the matrix algorithm, available at http://arxiv.org/abs/physics/0102055). 
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Appendix B: Divergence of the fonrth order dissipation source terms 

Again following Robertson , the general expression for a fifth order tensor of two-point 
type is given by 


Bniikl Bi 


rn ri Tj Tk ri 


f f j' 


^ ^ rk'f'i ^ ^ Thfi „ ^ r,- r, n ^ ^ r,- r,- ru „ , r„ ru n 

+ B2Sni -h B^Snj -h B^dnk -1“ B^6nl -h B^dij - 

'' rj^ rj^ rj^ rj^ rj^ ly^ ly '' ry^ ry^ ryt 


ry^ ryt ry^ 


f f y' rp ry^ rp 

Y' <7^ 7^, 'Y' 7^ 7^, lY' /Y' 7^ j /Y' 7^ 7^. 

+ + B.SU-J-- + + -Sii5, 


ry^ ry^ y^ ry^ ry^ ry^ ry^ ry' ryt 

Ti T]^ T • Ti 

-\- B\2^ni^jk H“ B\^5ni^jl B\^6ni5kl H“ B\^5'(ij6il^ “h B\Q6nj^il 

lyt y ry '' y '' y 

y . 'k'T' ' T' ' ^ k 

BljSnjSf^l “h Bi^S-fik^ij BigS-fikSil B2oSnkSjl “1“ B2lSfil6ij 

y y y y y 

y . y . y y y 

+ B22^nl^ik -1" B2'i5nl5jk -h B246ij6kl -h B2^5ik5ji -h B2eSil6jk —. 


rjT’ ip ip 


llOkl- 

Tk 


rp T' ■ T' ■ 

' n '2 ' j 


y y y 


(Bl) 


We consider tensors of the form {En^ijki) = {‘^UnEijki), 


i^^UnEijkl) {Skl E ^ku') “t” ^Ui^Uk i^^jl E E ^Ui^Ul i^^jk “t” ^jk'^ 

+AujAuk {eu + Sill) + AujAui {etk + Sik') + AukAui {sij + Eiji)]) , 


(B2) 


where the indices i, j, k and I may be interchanged, but n may not. Consequently, Bi = Ei, 
B 2 = ... = B^ = E 2 , Bq = ... = Bii = £'3, Bi 2 = ... = B 23 = E 4 ^ and £24 = £25 = B 2 Q = £ 5 . 
Thus we end up with hve scalar functions of r instead of three as we may not interchange 
all indices. That is, we may not use the result of Hill for the divergence of the fourth 
order structure functions. Next, we need to determine £j. We choose 


(Atti£4^o) = (B' 1 , 1111 ) = £1 -|- 4£2 -|- 6 E 3 + 12£4 -|- 3E^ 

(Arti£2,2) = (£1,1122) = -£'3-1- 2£4 -|- £5 

(Ami£o,4) = (£1,2222) = 3£5 

(Am2£3,i) = (£2,2111) = £2 + 3£4 

(Am2£i,3) = (£2,1222) = 3£4 

and solving for the £j gives 

£1 = (£1,1111) — 6 (£1,1122) + (£1,2222) — 4 (£2,2111) + 4 (£2,1222) 

£2 = (£2,2111) ~ (£2,1222) 

2 1 

£3 = (£1,1122) ~ - (£2,1222) — - (£1,2222) 

£4 = - (£2,1222) 

£5=2 (£1,2222) • 


(B3) 


(B4) 
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Taking the derivative of eq. fIBip using the relations eq. flA7j) then results in 


A 

dr„ 


' 7^ CjT' T' } 

2 ^ dEs 2 ^ dEi 2 ^ 

-E 2 H—^- 1 — E^ + 2 —-i ?4 

7^ CjT' 7^ Cyi^' T' 

r 'f'i'f'l r ^iVk ^ TiTj 

E5jk—Y E Sji —-—h 6ki—Y 

4 QE^ 2 \ 

--E 4 H—— 1 “ ~E^ 1 {Sij6ki + SikSji + Sii6jk) ■ 


Tkri , e- 


J.2 


+ ^ik — 5 “ + 


TjTk 


(B5) 


Therefore using eq. fIBSp . 


d (T^n,llll) 

d (7?!, 1111 ) 

dr 

dr 

d {En,1122) 

d (7?!, 1122 ) 

dr 

dr 

d {En, 2222 ) 

d (71^1,2222) 

dr 

dr 


{E 

{E 

{E 


1,1111 

1,1122 

1,2222 



(A' 2,2111 


) + -{E 

r 

4 

) + -{E 
r 


2,2111 

2,1222 


) 

)- 


), 


8 

3r 


(A' 2 , 1222 ) 


(B6) 


where 

(A'l^iiii) = (Atti£'4^o) = 6 (Att^ (^ii + £ 44 )) (B7) 

{E, ,1122) — {A.UiE 2 ^ 2 ) — ( Am4 [AM4 (^22 + ^22) E 4AM4M2 (£12 + £42) + ^^2 (£ll + £44)] ) 

(B8) 

{El ,2222) — {AuiEo^4) = 6 (AU4AM2 (£22 + £22)) (B 9 ) 

(-E'2,2111) = (Am2£'3,i) = (Am 2 [3AK2AK4 (£44 + £44) + 3Art4 (£42 + £42)]) (BIO) 

(-E'2,122) = {Au 2 Ei^^) — (^Au 2 [3AM4AM2 (£22 + ^22) E 3 Am 2 (£12 + ^12)]) • (Bll) 
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